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Применение Технологии композиционного 
преодоления неопределенности к решению 
задач с ограничениями 


В статье предложен метод решения задач линейного программирования, базирующийся на эволюционной 
парадигме. Исследован алгоритм его реализации, в основе которого находится полное пространство поиска 
возможных решений. Рассмотрены аспекты программной реализации метода. Выполнена экспериментальная 
верификация и сравнительный анализ с результатами применения генетических алгоритмов. 


Введение 


Функционирование сложных производственных, экономических, технических 
и социальных систем характеризуется набором определенных факторов, параметров 
и оценивается множеством критериев. Как факторы, так и критерии могут иметь 
различную природу, но в подавляющем большинстве случаев (в практических задачах), 
количество критериев ограничено, а на значения факторов или их комбинации наложены 
определенные ограничения. Математически формализованные задачи оптимизации 
процессов функционирования указанных систем составляют класс задач с ограни- 
чениями (Сопзгаше Зан$Расйоп Ртоетл) [1]. 

К классу задач с ограничениями относятся и задачи линейного программирования, 
широко применяемые в экономико-математическом моделировании. Использование 
таких моделей в экономике позволяет выделить и формально описать наиболее 
существенные связи между экономическими объектами, используя высокий уровень 
абстракции в силу сложности изучаемых процессов и явлений, индуктивным путем 
получить новые знания об объекте исследования. Первоочередными проблемами, для 
решения которых целесообразно применять экономико-математическое моделирование 
и методы оптимизации, являются: определение номенклатуры продукции (видов услуг), 
а также оптимизация объектов производства на определенный перспективный период; 
распределение имеющихся материальных и финансовых ресурсов по видам деятель- 
ности; определение цены, которая будет обеспечивать оптимальный уровень прибыли; 
накопление собственных средств для развития производственной деятельности, опре- 
деления потребностей в кредитах, которые могут привлекаться; определение влияния 
изменений стоимостных показателей на экономическую эффективность предприятия и т.д. 
Однако решения таких задач является сложным нетривиальным процессом, который 
требует разработки эффективных методов поиска решений с использованием компью- 
терной техники, новых идей, принципов, моделей и методов. 

В последние годы в области эволюционных вычислений новые методы появляются 
из идеи комбинированного использования технологий Зой Сотрийпо [2]. Одним из 
таких методов является технология композиционного преодоления неопределенности [2], 
которая позволяет осуществлять целенаправленную оптимизацию, учитывающую 
субъективные предпочтения исследователя и интегрирующая в себе элементы 
теории вероятностей, теории нечетких множеств и эволюционных стратегий. 
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Целью настоящей работы является адаптация технологии композиционного 
преодоления неопределенности к решению задач линейного программирования на 
основе использования штрафных функций, определение особенностей ее применения, 
преимуществ и недостатков, а также выполнение сравнительного анализа результатов, 
полученных с использованием эволюционных методов. 


Постановка задачи 
Задачи С ограничениями заключаются в поиске: 


у= (Хор, хЕХ" (1) 
при ограничениях 
Мм 


с= (*.,5. СЕС”", с= С: (2) 


т=1 

где Х = (еде ] — множество переменных, принадлежащих пространству 
Е (х) — целевая функция; г, < Х” -— произвольное т -арное отношение на Хх; 
5, :Х" > Х” — функция для проектирования вектора Х’ = (х,,...х, ) Е Х" на некоторые 
из т его компонент [1]. 

Классифицируют задачи с ограничениями, как правило, по свойствам входной 
информации (детерминированные, стохастические); по характеру описываемых явлений 
(динамические, статистические); типу переменных (зависимые, независимые; дискрет- 
ные, непрерывные; детерминированные, случайные); типу ограничений (унарные, 


бинарные, высших порядков; линейные, нелинейные; абсолютные, приоритетные). 
Задачи линейного программирования формулируются так: 


у= Л (Х) > орь хЕХ" (3) 

при ограничениях 
&(Х)=0,1=0,...р, р>0, (4) 
и (Хх) йе р+ь,..т, т-р>0, (5) 


где хЕЗПР. Множество 5 < Х" определяет пространство поиска, а множество 
Е < Х" - допустимое пространство поиска, добавленное ограничениями. Пространство 


поиска 5 является и-мерным прямоугольником в Х" (область определения перемен- 
ных определена верхними и нижними границами): 


КО <х, <и(1, 1<1< и. 


Основные направления решения задачи 


Поскольку необходимо найти глобальные оптимумы целевой функции в много- 
мерном пространстве независимых переменных, целесообразно применить эволюци- 
онные методы. Их использование, в отличие от традиционных методов, не зависит от 
выбора начальной точки и не требует выполнения дополнительных условий на 
характеристики целевой функции. Выбор метода эволюционного моделирования 
является прерогативой исследователя. 

Для обработки ограничений с целью изъятия неперспективных решений в 
генетических алгоритмах применяются штрафные функции: 


Х Хх}, если ХЕЁЕП5, 
еуа1(Х) = } 


Г(х)+ репапу (х), в противном случае, 
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где штраф репайу(Х ) равен нулю, если ограничение выполняется, в противном 
случае штраф имеет определенное значение. Существует несколько подобных 
методов, в большинстве из них множество функций Л, (1 = т) используется для 


генерации штрафов; функция /, оценивает штраф, который будет наложен на /-е 


ограничение, по формуле: 


тах, 5, [99 якщо 1< ]<4 


7) и (х] якщо а+1< } < т. 


Конечно, методы отличаются некоторыми важными деталями, в частности, построением 
штрафных функций и изъятием неперспективных решений. Рассмотрим два из них. 

Первый метод, предложенный А. Нота!#г, 5.Н.-У. Гав Х. 01 [3], предусматривает, 
что для каждого ограничения формируется множество интервалов, которые определяют 
соответствующие штрафные коэффициенты. Метод представлен таким алгоритмом: 

Шаг 1. Для каждого ограничения сформировать несколько (1 ) штрафных уровней. 

Шаг 2. Для каждого штрафного уровня и каждого ограничения сгенерировать 
штрафные коэффициенты К, ( =1.2...,1,1=12,..., т) ‚ причем наивысшему штрафному 
уровню должно соответствовать наибольшее значение коэффициента. 

Шаг 3. Определить выборочную совокупность возможных решений. 

Шаг 4. Вычислить значение функции, оптимумы которой ищут, в точках 
выборочной совокупности по формуле: 


еха(Х) = Г(х)+ У, ВЛ, (х) 


Недостаток метода заключается в том, то он зависит от большего числа 
параметров. Так, для т ограничений метод требует: т параметров для генерации 
числа интервалов для каждого ограничения — границ интервалов, / параметров для 
каждого ограничения — штрафных коэффициентов К,. Таким образом, в общем случае 


метод требует т(21 + 1) параметров для управления т ограничениями. 


Авторы второго метода /.А. Лютез и С.В. НоиК [4], аналогично предыдущему 
методу, предлагают использовать динамические штрафы. Индивиды вычисляются 
(на итерации 1) с использованием формулы: 


еха(Х | = Г(х)+ (сх!) уе А () 


где С, х и ДВ константы. Целесообразно принять такие значения этих параметров: 
С =0,5, а = В = 2. Этот метод содержит намного меньше параметров, чем предыдущий. 
Воспользуемся приведенными штрафными функциями в технологии компози- 
ционного преодоления неопределенности [2], которая определена таким алгоритмом 
(процедура Еуотах): 
Шаг 1. 1=1. 
Шаг 2. Определить совокупность возможных решений (=) и определить 


процентное соотношение количества индивидов на следующем шаге. 
Шаг 3. Пока условия завершения работы алгоритма не выполнены { 
Шаг3.1. 1=1+1. 


Шаг 3.2. Вычислить значение функции(=). 
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Шаг 3.3. Построить функцию принадлежности /м. 
Шаг 3.4. Определить множество точек {х,}, которые принадлежат множеству 


среза и для которых выполняется неравенство , _ х, НО". 
Шаг 3.5. Генерировать последовательности {2 } для каждой точки с 


множества {х, }. 
Шаг 3.6. Мутация а 
} 


Приведенный метод позволяет осуществлять целенаправленную оптимиза- 
цию, учитывая цели и задачи исследования, поскольку в нем в значительной степени 
преодолены такие недостатки технологий Зой Сотрийпе, как: зависимость точности 
решений от дискретности их представления и заданной точности будущего результата 
(как в генетических алгоритмах); использование исключительно равномерного и (или) 
нормального распределения, что значительно увеличивает время решения задачи 
(как в эволюционных стратегиях); зависимость сходимости от качества популяции. 


Результаты экспериментов 


Рассмотрим пример. Фабрика производит изделия четырех наименований: печенье, 
пряники, сухари и бублики. Для изготовления продукции используется сахар, расти- 
тельное масло, пшеничная мука, маргарин, изюм. 

Интенсивность потребления каждого ингредиента на тонну изделия соответству- 
ющего наименования и прибыль от реализации единицы веса заданы в табл. 1. 
Необходимо определить оптимальный объем производства продукции каждого наимено- 
вания (на фабрике есть 20 т сахара, 24 т растительного масла, 30 т муки, 0,8 т 
маргарина, 0,6 т изюма), при котором прибыль будет максимальной. 


Таблица 1 — Интенсивность потребления ингредиентов и прибыль 


Ингредиенты 

Сахар Раститель Мука Марга- И Прибыль, 

И ное масло рин тыс. грн 
зделие 

Печенье 0,15 0,2 0,6 0,03 0,02 1 
Пряники 0,4 0,1 0,8 — — 1,2 
Сухари 0,3 0:2 0,5 —- 0,05 0,8 
Бублики 0,2 0,2 0,4 — — 1 


Рыночный спрос на печенье, пряники, сухари и бублики составляет 5, 6, Зи 4т 
соответственно. 


Математическая модель задачи. Пусть х, > 0,1= 13 — количество изделий 
1-го наименования, тогда экономико-математическая модель задачи выглядит так: 
0,15х, + 0,4х ,+0,3х, + 0,2х, < 20, 
0,2х, + 0,1х,-+0,2х, + 0,2х, < 24, 
0,6х, +0,8х,-+0,5х, + 0,4х, < 30, 
0, 0Зх, < 0,8, 
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0,02, + 0,05х, < 0,6, 
Оо, 
О <х, < 6, 
= 
О<х, <4, 
Д=х +Ъ2х, +0,8х, + х, >> шах. 

Для исследования, целью которого было выяснение эффективности технологии 
композиционного преодоления неопределенности при решении задач линейного 
программирования, выполнено фиксированное число запусков алгоритма. При этом 
начальная популяция состояла из 50 индивидов, максимальное число возможных 


итераций — 300, вероятность мутации — 1%, среднеквадратическое отклонение - 1, 
точность решения — 10“. 


Таблица 2 — Результаты экспериментов 


Генетический Композиционный 
Название критерия алгоритм метод 
Метод | | Метод2 | Метод1 | Метод 2 
. Ь | 16,767 | 91618 18,6 18,6 
Критерий ш| 9,5754 | 12871 | 12,655 | 12,428 
максимума | 13,157 | 17869 | 15,128 | 15,395 
функции { 16,6 35 124 4,46 6,096 
№ з ь 18,6 18,6 
я а. т . . 12,63 12,883 
Е. нае 14/769 | 15,691 
а { 107,109 | 46,561 
Критерий Ь | 21,606 | 18,059 18,6 18,6 
максимального ш | 12,820 9,17 12,785 12,427 
количества \ | 16,968 13,625 15,664 15,07 
итераций { 88,45 86,53 6,46 6,26 
Ь | 21,679 21,646 18,6 18,6 
ш| 12,707 | 12,758 | 12,509 | 12,594 
ИРИ | 17079 17879 | 14/2 | 15,603 
| 50,21 79,42 3,42 2,43 


Результаты экспериментов для каждого метода (1, 2) поданы в табл. 2, где 
показано наилучшее (Ъ), среднее (т) и наихудшее (\) значение целевой функции, а 
также время расчета (1). Эксперименты проводились с использованием трех типов 
критериев. Согласно первому критерию итерации прекращались, когда максимальные 
значения целевой функции на соседних итерациях незначительно отличались. По второму 
критерию сравнивались средние значения целевой функции на разных итерациях. 
Вычисления прекращались после определенного количества итераций — по третьему 
критерию. Использовался также интегральный критерий: вычисления прекращались, 
если выполнялся хотя бы один критерий из предложенных выше. 

Выводы, которые можно сделать по результатам моделирования, являются неодно- 
значными. В частности, согласно критерию среднего значения функции, используя 
генетический алгоритм, не удалось получить результат за приемлемое время. Лучшие 
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результаты были получены в случае применения критерия максимального количества 
итераций. Причиной этого является полиэкстремальный характер целевой функции 
и, как следствие, частое попадание в локальные экстремумы. 


Выводы 


Предложенный метод позволяет решать задачи линейного программирования, 
учитывая особенности ограничений. Еще одним его преимуществом является время 
решения задачи. Как следует из табл. 2, значение времени, необходимое для расчета 
разработанным методом, на порядок меньше времени работы генетического 
алгоритма. Результаты работы двух методов являются сравнимыми. 

Следует заметить, что технология композиционного преодоления неопределен- 
ности является параметрическим методом. Его использование для решения указанной 
задачи осуществлено без оптимизации параметров. Применение метода для решения 
тестовых оптимизационных задач показало, что оптимизация параметров метода 
приводит к увеличению точности результата до 30% и скорости сходимости до 24%. 
Поэтому параметрическая оптимизация предложенного метода является еще одним 
перспективным способом повышения эффективности процессов решения задач с 
ограничениями, тем более что для эволюционных методов эта проблема актуальна. 
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Аспекти застосування методу композищИного подолання невизначеност! в задачах з обмеженнями 

У статт! запропонований метод розв’язання задач лиЙного програмування на основ! еволющйнот 
парадигми. Досл1джено алгоритм його реал1заци на основ! повного простору пошуку можливих ришень. 
Розглянуто аспекти програмно! реалзацй методу. Виконана експериментальна вериф\кащя та наведен! 
результати пор1вняльного анал1зу з генетичними алгоритмами. 


О.Т. Уезотоъа, Г.Е. Зпури 

Аррйсабоп ое Опсенашеу Сотрозвоп Оуегсотшо ТесВпоосу ® Зауше РгоШет$ уу Соп$гат6 

№1 5 рарег Фе рго Мет 4ес151оп тейо4 юг Ппеаг ргоэтатите оп Фе Ъаз15 оЁ еуоаНопагу рага@ ет 1$ 
оегед. № геайтаЯюоп а]еогит 1$ шуезИзайопа! оп йе Ба$15 оЁ сотр/ее зрасе зеагсь о{ розз1е 4ес110п5. 
ТБе азрес5 оЁ тефо4 зойу\аге геаПхайоп аге сопз1еге4. Ехрегипета| уепйсаноп ап4 сотрагайуе апа[уз1$ 
фе гези[5 аге ехесше4. 


Статья поступила в редакцию 22.06.2011. 
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